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Abstract. We present global fits of the constrained Minimal Supersymmetric Standard 
Model (cMSSM) and the Non-Universal Higgs Model (NUHM), including the most recent 
CMS constraint on the Higgs boson mass, 5.8 fb _1 integrated luminosity null Supersymmetry 
searches by ATLAS, the new LHCb measurement of BR(B S — > /i + '//") and the 7-year WMAP 
dark matter relic abundance determination. We include the latest dark matter constraints 
from the XENON100 experiment, marginalising over astrophysical and particle physics uncer- 
tainties. We present Bayesian posterior and profile likelihood maps of the highest resolution 
available today, obtained from up to 350M points. We find that the new constraint on the 
Higgs boson mass has a dramatic impact, ruling out large regions of previously favoured 
cMSSM and NUHM parameter space. In the cMSSM, light sparticles and predominantly 
gaugino-like dark matter with a mass of a few hundred GeV are favoured. The NUHM ex- 
hibits a strong preference for heavier sparticle masses and a Higgsino-like neutralino with a 
mass of 1 TeV. The future ton-scale XENON1T direct detection experiment will probe large 
portions of the currently favoured cMSSM and NUHM parameter space. The LHC operating 
at 14 TeV collision energy will explore the favoured regions in the cMSSM, while most of 
the regions favoured in the NUHM will remain inaccessible. Our best-fit points achieve a 
satisfactory quality-of-fit, with p-values ranging from 0.21 to 0.35, so that none of the two 
models studied can be presently excluded at any meaningful significance level. 
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1 Introduction 

The Large Hadron Collider (LHC) is providing valuable information on possible extensions 
of the Standard Model of particle physics. In particular, LHC data are probing some of the 
most popular models of Supersymmetry (SUSY). At the same time, dark matter experiments 
searching for direct evidence of dark matter particles scattering off atomic nuclei are providing 
stringent constraints on the parameter space of supersymmetric models, that are highly 
complementary to the LHC results (see e.g. Refs. [1] and references therein). 

We have recently presented global fits of one of the most popular SUSY models, the 
so-called constrained Minimal Supersymmetric Standard Model (cMSSM, see e.g. [2, 3]), 
including LHC data with 1 fb _1 of integrated luminosity and XENON100 direct detection 
data presented in 2011 [4, 5]. Other global fits studies of the cMSSM can for example be 
found in Ref. [6-13]. In this paper, we extend our analysis to the Non-Universal Higgs Model 
(NUHM, see e.g. [14-16]), a less constrained model of SUSY with less restrictive boundary 
conditions applied at the Grand Unification scale than in the cMSSM. 

In the following we discuss the impact of new LHC data, including the discovery of 
the Higgs boson, and the latest XENON100 data on both the cMSSM and the NUHM. 
More specifically, we implement in our cMSSM global fits the latest exclusion limits in the 
(mo,m 1( / 2 ) plane obtained by the ATLAS collaboration for squarks and gluinos in final states 
that contain missing Et, jets and leptons in 5.8 fb _1 integrated luminosity of data at 
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y/s = 8 TeV collision energy [17]. For the NUHM we additionally include the exclusion limit 
in the (m^tan^) plane from a CMS search for the decay of neutral Higgs bosons into final 
states containing two muons and missing Et, based on 4.5 fb -1 integrated luminosity of data 
collected at yfs = 7 TeV collision energy [18]. 

Furthermore, we include the most recent experimental constraint from the CMS collab- 
oration on the mass of the lightest Higgs boson rrih = 125.8 ± 0.6 GeV. This constraint is 
derived from a combination of 5.1 fb" 1 yfs = 7 TeV data and 12.2 fb" 1 y/s = 8 TeV data [19]. 
For the first time, we include the new measurement of BR(B S — > = (3. 2^^) x 1CP 9 

from the LHCb collaboration, derived from 1 fb -1 of data at \fs = 7 TeV collision energy 
and 1.1 fb -1 of data at y/s = 8 TeV collision energy [20]. Finally, we also include the most 
recent dark matter results from the XENON 100 direct detection experiment, obtained from 
224.6 live days of data and 34 kg fiducial volume, collected between February 2011 and March 
2012 [21]. 

The main aim of this paper is to study the impact of new experimental results from LHC 
SUSY and Higgs searches, from the LHCb collaboration, and from direct detection searches 
for dark matter with 225 live days of XENON100 data on the cMSSM and the NUHM. We 
determine the most favoured regions of the cMSSM and NUHM parameter space in light of 
these experimental constraints, from both the Bayesian and the profile likelihood statistical 
perspective. We discuss the viability of the best-fit points for both models, and evaluate the 
extent to which these models are disfavoured by present-day experimental constraints. 

The paper is organized as follows. In section 2 we present the theoretical and statistical 
framework, and provide details about the supersymmetric models we study, the priors used for 
the model and nuisance parameters, and the implementation of the experimental constraints. 
In section 3 we present the favoured regions and best-fit point for the cMSSM and discuss 
the impact of different experimental constraints, in particular of the g — 2 measurement, as 
well as the prospects for discovery with future direct detection experiments and the LHC 
operating at 14 TeV collision energy. In Section 4 we repeat this analysis for the NUHM, and 
discuss the phenomenological differences between the cMSSM and the NUHM. We present 
our conclusions in section 5. 

2 Theoretical and statistical framework 

In this paper we consider two different SUSY models, the cMSSM and the NUHM. The 
cMSSM is a model that has been widely studied in the past. It employs strong universality 
conditions at the Grand Unification scale, as a result it can be described by only five free 
parameters: the universal scalar and gaugino mass parameters mo and JTfx/2, the universal 
trilinear coupling Aq, the ratio of the two Higgs vacuum expectation values tan f3, and the sign 
of the Higgs/Higgsino mass parameter sgn(/i). In the following analysis we fix sgn(//) = +1, 
which is favoured by the measured value of the magnetic moment of the muon [22]. The 
universality assumptions in the cMSSM are motivated by the natural link between SUSY 
and Grand Unification Theories; due to its low number of parameters the cMSSM has been 
an extremely popular model to study in the past. 

A more general model with less restrictive boundary conditions applied at the Grand 
Unification scale is the NUHM. Specifically, in the NUHM the soft masses of the two Higgs 
doublets m,H u and mff d are free parameters, independent of the sfermion masses. This is 
a reasonable assumption to make, since the Higgs and the matter fields are described by 
different supermultiplets, so that there is no strong motivation to assume unified Higgs and 
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sfermion masses. The NUHM is described by six free parameters: the four cMSSM parameters 
mo, Tii/2, A) j tan/3, the (continuous) Higgs/Higgsino mass parameter \i and the mass of the 
pseudoscalar Higgs, tua (instead of the last two parameters one could use mn u and rriH d , 
as has been done in previous works, e.g. [23]). The larger number of free parameters in the 
NUHM compared to the cMSSM leads to a richer phenomenology. 

Although we use highly efficient Bayesian methods to explore the favoured regions of the 
cMSSM and NUHM parameter spaces given all available constraints, we present our results 
both in Bayesian and in frequentist terms. Our approach starts from Bayes' theorem [24] 

,(e,D) - („ 

where D are the data and are the model parameters of interest. Bayes' theorem states that 
the posterior probability distribution function (pdf) p(0|D) for the parameters is obtained 
from the likelihood function p(D|0) = £(0) and the prior pdf (or "prior" for short) p(&). 
In this paper we are primarily interested in parameter inference, therefore the Bayesian 
evidence p(D) merely acts as a normalisation constant, and will not be considered further in 
the following analysis. 

In order to study the constraints on a single parameter of interest 6i, one can consider 
either the one-dimensional marginal posterior pdf (Bayesian), or the one-dimensional profile 
likelihood (frequentist). The marginal pdf is obtained from the full posterior distribution by 
integrating (marginalising) over the unwanted parameters in the n-dimensional parameter 
space: 

p(0i\D) = f P {®\D)de 1 ...d6 i . l dd l+1 ...d6 n . (2.2) 

The frequentist profile likelihood function for (9j, instead, is found by maximising over 
the parameters that are not of interest: 

£(0i) = max £(0). (2.3) 

01, — ,8i—l)8i+l, — ,6 n 

The extension of these concepts to more than one parameter is straightforward. The profile 
likelihood function and the marginal posterior pdf are two different statistical quantities 
that may lead to different conclusions about the parameter space of interest. The marginal 
posterior pdf integrates over hidden parameter directions and therefore correctly accounts for 
volume effects; it peaks at the region of highest posterior mass. The profile likelihood function 
peaks at the region of highest likelihood. It is oblivious to volume effects, but is an excellent 
quantity to find small regions of high likelihood in parameter space. Especially when studying 
complicated parameter spaces of high dimensionality, such as SUSY parameter spaces, these 
two quantities will usually not lead to the same conclusions, and the maximum of information 
about the model parameter space is obtained by studying both of these quantities [25, 26]. 
Therefore, in the following we present results for both the marginalised Bayesian posterior 
pdf and the profile likelihood function. 

2.1 Model and nuisance parameters 

As can be seen from Eq. (2.1), the posterior pdf is dominated by the likelihood whenever the 
prior is flat in the parameter space region where the likelihood has support. However, this 
is not generally the case, and therefore the choice of prior distribution can have a significant 
impact on the resulting posterior. One approach is the adoption of an informative prior, that 
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cMSSM Parameters 




Flat priors Log priors 


cMSSM parameters 


m [GeV] 
m x / 2 [GeV] 
A [GeV] 


(50.0, 4000.0) (10 i 7 , 10 3 - e ) 
(50.0, 4000.0) (10 17 , 10 3 - 6 ) 
(-4000.0, 4000.0) 


tan /3 


(2.0, 65.0) 


NUHM parameters as above, and additionally: 


H [GeV] 
m A [GeV] 


(-2000.0, 2000.0) 
(50.0, 4000.0) (10 1 - 7 , 10 3 - 6 ) 



Table 1. cMSSM and NUHM parameters and the range of their values explored by the scan. Flat priors are 
uniform in the masses; log priors are uniform in the logarithm of the masses. 

for example encapsulates in a natural way the theoretical prejudice that finely-tuned regions 
ought to be penalised (in an Occam's razor sense, see e.g. [24]). In this work, we compare 
the posterior distributions for two different choices of priors, in order to evaluate the prior 
dependence of the posterior. Such a prior-dependence has been found to be commonplace in 
high-dimensional parameter spaces with complex, multi-modal likelihoods which are typical 
of SUSY phenomenology (see e.g. [4, 5, 27]). If the posterior distribution displays some 
residual dependence on the choice of prior distribution, the resulting constraints on the 
observables have to be interpreted with care [26, 27]. 

In order to assess to what extent the posterior distribution is influenced by the choice 
of prior, we repeat each of our scans for two different choices of (non-informative) prior 
distributions. "Flat" priors are uniform on the mass parameters of the cMSSM (mo, m^) 
and NUHM (mo, mjrt, m,^), "log" priors are uniform in the log of the mass parameters. 
Both choices of priors are uniform on Aq, tan/3 (cMSSM and NUHM) and \i (NUHM). Prior 
ranges for the cMSSM and the NUHM are summarised in Table 1. 

In addition to the above model parameters, we include several nuisance parameters in 
the scans. Residual uncertainties on measurements of certain Standard Model parameters 
have been shown to have an important impact on the results of SUSY studies [28]. To 
correctly account for this effect, we include four SM parameters as nuisance parameters 
in our scan (the top mass, the bottom mass, the electroweak coupling constant and the 
strong coupling constant). When including direct detection constraints in the analysis one 
additionally needs to take into account uncertainties in astrophysics and the nuclear physics 
quantities. Here, we adopt the same strategy as presented in Ref. [5], and include three 
additional hadronic nuisance parameters (the hadronic matrix elements /r u , f?d an d Jts 
that parameterise the contribution of the light quarks to the proton composition) and four 
astrophysical nuisance parameters (the local dark matter density p\ oc , and three quantities 
parameterising the WIMP velocity distribution: the local circular velocity v\ ST , the escape 
velocity -y esc and the velocity dispersion Vd) in the analysis, see [5] for full details. Nuisance 
parameters are well constrained by experimental data, and therefore we adopt informative 
Gaussian priors for these quantities. The mean and standard deviation of the Gaussian priors 
are chosen to reflect up-to-date experimental constraints, and are given in Table 2. 
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Gaussian prior 


Range scanned 


Ref. 


SM nuisance parameters 


Mt [GeV] 


173.2 ±0.9 


(170.5, 175.9) 


[29] 


mb(mb) M 


[GeV] 4.20 ± 0.07 


(3.99, 4.41) 


[30] 


[a em {M z ) MS ]- x 127.955 ± 0.030 


(127.865, 128.045) 


[30] 


a s {M z ) MS 


0.1176 ±0.0020 


(0.1116, 0.1236) 


[31] 


Astrophysical nuisance parameters 


pioc [GeV/cm 3 ] 0.4 ±0.1 


(0.1, 0.7) 


[32] 


v hr [km/s] 


230.0 ±30.0 


(140.0, 320.0) 


[32] 


u csc [km/s] 


544.0 ± 33.0 


(445.0, 643.0) 


[32] 


v d [km/s] 


282.0 ±37.0 


(171.0, 393.0) 


[32] 


Hadronic nuisance parameters 


fTu 


0.02698 ± 0.00395 


(0.015, 0.039) 


[33] 




0.03906 ± 0.00513 


(0.023, 0.055) 


[33] 


fTs 


0.363 ±0.119 


(0.0006, 0.72) 


[33] 



Table 2. List of nuisance parameters included in the scans of the cMSSM and NUHM parameter spaces. 
The mean value and standard deviation adopted for the Gaussian prior on each parameter is shown, as well 
as the range of values explored by the scan. 



2.2 Experimental constraints and the likelihood function 

The likelihood function entering in Eq. (2.1) is composed of several different parts, corre- 
sponding to the different experimental constraints that are applied in our global fits analysis: 

ln£ = ln£ LH c + hi£ W MAP + ln£ E w + ln£ B (D) + ln£ 9 -2 + m£xcioo- (2.4) 

The full list of experimental constraints included in the likelihood function is given in 
Table 3. 

• In £lhc • The LHC likelihood implements recent results from SUSY null searches from 
ATLAS and CMS. Exclusion limits in the (?no> m i/2) plane are based on a search by the 
ATLAS collaboration for squarks and gluinos in final states that contain missing , 
jets and leptons in 5.8 fb ~ l integrated luminosity of data at sfs = 8 TeV collision 
energy [17]. While this exclusion limits was obtained in the cMSSM framework for 
fixed values of tan/3 = 10 and Aq = 0, the result is fairly insensitive to tan/3 and Aq, 
so that we can use this limit even when varying these quantities. Furthermore, we 
can also apply this limit to the NUHM, since the signal is dominated by the strong 
cross-sections of the two first generation of squarks and gluino production, which do not 
depend significantly on other parameters than ttiq and ra\ji . The LHC exclusion limit 
is included in the likelihood function by defining the likelihood of samples corresponding 
to masses below the limit to be zero. For the NUHM we additionally include the 
exclusion limit in the (mA,tan/3) plane from a CMS search for the decay of neutral 
Higgs bosons into final states containing two muons and missing Er, based on 4.5 fb -1 
integrated luminosity of data collected at y/s = 7 TeV collision energy [18]. 

1 Additionally, as we will show below, the LHC exclusion limit has essentially no impact on the NUHM 
results, so that small variations of this exclusion limit will not change our conclusions. 
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We furthermore include the most recent experimental constraint from the CMS collab- 
oration on the mass of the lightest Higgs boson = 125.8 ± 0.6 GeV. This constraint 
is derived from a combination of 5.1 fb _1 = 7 TeV data and 12.2 fb _1 y/s = 8 TeV 
data [19]. The statistical significance of the signal is 6.8cr. We use a Gaussian likelihood 
and we add in quadrature a theoretical error of 2 GeV to the experimental error. We 
do not impose the experimental constraint on the Higgs production cross-section in this 
analysis, since all of our samples fall within a very narrow range ct/ctsm = [0.95, 1.00], 
which is in good agreement with the CMS constraint ct/ctsm = 0.88 ± 0.21 reported in 
[19]. 

Finally, we include the new LHCb constraint on BR(B S — > = '§) x 10 -9 , 

derived from a combined analysis of 1 fb _1 data at \fs = 7 TeV collision energy and 
1.1 fb" 1 data at \/s = 8 TeV collision energy [20]. We implement this constraint as a 
Gaussian distribution with a conservative experimental error of a = 1.5 x 10~ 9 , and a 
10% theoretical error. 

• 1ii£\vmap- The WMAP measurement of the dark matter relic abundance is included as 
a Gaussian in ln£\yMAP- We use the WMAP 7-year value £l x h 2 = 0.1109 ± 0.0056 and 
we add a fixed 10% theoretical uncertainty in quadrature. We assume that neutralinos 
make up all of the dark matter in the universe. 

• 1ii£ew implements precision tests of the electroweak sector. The electroweak precision 
observables M\y and sin 2 9 e ff are included with a Gaussian likelihood. 

• ln£g(rj). The B and D physics observables are included as a Gaussian likelihood. The 
full list of B and D physics observables included in our analysis is shown in Table 3. 

• ln£ 9 _2- The measured anomalous magnetic moment of the muon, included as a Gaus- 
sian datum, provides important information about supersymmetric parameter spaces, 
since it can be experimentally measured to very good precision. By comparing the the- 
oretical value of this quantity favoured in the Standard Model with the experimental 
result [35] the supersymmetric contribution 5a^ USY to this quantity can be constrained. 
The experimental measurement of the muon anomalous magnetic moment is discrepant 
with the SM prediction by 5a^ USY = (28.7 ± 8.0) x 10~ 9 [35]. This analysis, based on 
e + e~ data, leads to a 3.6<r discrepancy between the experimental result and the theoret- 
ical SM value. However, the corresponding analysis based on r data leads to a smaller 
discrepancy of 2.4<r [35]. The significance of this discrepancy has to be interpreted with 
care, since the calculation of the theoretical value of the muon anomalous magnetic 
moment in the Standard Model is subject to important theoretical uncertainties, aris- 
ing mostly in the computation of the hadronic loop contributions. In past studies of 
the cMSSM [4] and NUHM [23] it was found that the 5a s ^ SY constraint plays a domi- 
nant role in driving the global fits results for both models. Therefore, in the following 
we will present results for both an analysis including and excluding the experimental 
constraint on <5a^ USY , in order to evaluate the dependence of our conclusions on this 
constraint. 

• ln£xcioo- We also include the most recent dark matter results from the XENON100 
direct detection experiment, obtained from 224.6 live days of data and 34 kg fiducial 
volume, collected between February 2011 and March 2012 [21]. XENON100 currently 
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places the tightest direct detection constraints on the WIMP properties. The collabora- 
tion reported the detection of two candidate WIMP scattering events in the pre-defined 
signal region, with an expected background of b = 1.0 ±0.2 events. The detected events 
are compatible with the background, so that new exclusion limits could be derived in 
the (m^o , a~l ) plane. 

We use an updated version of the approximate likelihood function from our previous 
analyses [4, 5], for a detailed description we refer the reader to Ref. [5]. In addi- 
tion to changes in the fiducial mass and exposure time, and the number of detected 
and background events, updates to the likelihood function include the reduction of the 
lower energy threshold for the analysis to 3 photoelectron events (PE) and an update 
to the response to 122 keV gamma-rays from calibration measurements to L y = 2.28 
PE/keVee, in accordance with the values reported in Ref. [21]. We make the simplify- 
ing assumption of an energy-independent acceptance of data quality cuts, and adjust 
the acceptance-corrected exposure to accurately reproduce the exclusion limit in the 
(m^,o, cj?o_ p ) plane reported in Ref. [21] in the mass range of interest. 

2.3 Scanning methodology 

To map out the posterior pdf and the profile likelihood we use the SuperBayeS-v2 . pack- 
age, an evolution of the publicly available SuperBayeS-vl . 5 [5, 27, 28, 42], which has been 
developed for this work. SuperBayeS-v2 . will shortly be released to the public. 

This latest version of SuperBayeS is interfaced with SoftSUSY 3.1.7 as SUSY spectrum 
calculator, MicrOMEGAs 2.4 [43, 44] to compute the abundance of DM, DarkSUSY 5.0.5 [45, 
46] for the computation of &^>_ p an d ^oL^ Superlso 3.0 [47, 48] to compute <5a^ USY and 

B(D) physics observables, SusyBSG 1.5 for the determination of BR(B — > Xg'y) [49, 50] 
and FeynHiggs 1.9 [51] to compute the Higgs production cross-sections. The infrastructure 
of the code has also been revised, to allow for a larger choice of super symmetric models at 
compilation time, including (but not limited to) the cMSSM and the NUHM studied here. 

The SuperBayeS-v2 . package uses the publicly available MultiNest v2.18 [52, 53] 
nested sampling algorithm to explore the cMSSM and NUHM model parameter space. Multi- 
Nest is an extremely efficient scanning algorithm that can reduce the number of likelihood 
evaluations required for an accurate mapping of the posterior pdf by up to two orders of 
magnitude with respect to conventional MCMC methods. This Bayesian algorithm, origi- 
nally designed to compute the model likelihood and to accurately map out the posterior, 
is also able to reliably evaluate the profile likelihood, given appropriate MultiNest settings, 
as demonstrated in [25]. We use the settings recommended in Ref. [25] (number of live 
points nlive = 20,000, tolerance tol = 10 -4 ), tuned to obtain an accurate map of the profile 
likelihood function. 

Our cMSSM posterior results are based on approximately 220M (128M) likelihood eval- 
uations for log (flat) priors. cMSSM posterior results excluding the constraint on <5a^ USY are 
derived from 199M (124M) likelihood evaluations for log (flat) priors. For the NUHM, the 
chains for log (flat) priors were generated from 132M (73M) likelihood evaluations. NUHM 
results excluding the Ja^ USY constraint are based on 95M (52M) likelihood evaluations. The 
profile likelihood function, which is in principle prior-independent, is derived from combined 
chains of the log and flat prior scans. In order to achieve a higher resolution even in the 
tail of the profile likelihood, we save the value and coordinates of all likelihood evaluations 
for the profile likelihood analysis. This procedure thus includes a large number of samples 
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Observable 


Mean value 


Uncertainties 


Ref. 






a (exper.) 


r (theor.) 




M w [GeV] 


80.399 


0.023 


0.015 


[34] 


sin 2 9 e f f 


0.23153 


0.00016 


0.00015 


[34] 


5a s v SY x 10 10 


28.7 


8.0 


2.0 


[35] 


BR{B ->• X s j) x 10 4 


3.55 


0.26 


0.30 


[36] 


Ram Bs 

1 i III 1 J ,1 }TV J 

BR(B u ^tu)sm 

A _ x 10 2 

BR(B-^Dru) 1n2 
RRl R^DeiA u 


1.04 
1.63 
3.1 

41.6 
0.999 


0.11 

0.54 

2.3 

12.8 

0.007 


- 
- 
- 

3.5 

_ 


[37] 
[36] 
[38] 
[39] 
[40] 


BR(D S tv) x 10 2 


5.38 


0.32 


0.2 


[36] 


5#(L> S //i/) x 10 3 


5.81 


0.43 


0.2 


[36] 


BR(D iiv) x 10 4 


3.82 


0.33 


0.2 


[36] 


n x h 2 


0.1109 


0.0056 


0.012 


[41] 


m h [GeV] 


125.8 


0.6 


2.0 


[19] 


BR{B S -)■ jU+aO 


3.2 x 10~ 9 


1.5 x 10" £ 


10% 


[20] 




Limit (95% CL) 


r (theor.) 


Ref. 


Sparticle masses 


As 


in table 4 of Ref. [42]. 




mo,m 1/2 


ATLAS, ^fs 


= 8 TeV, 5.8 fb 


_1 2012 limits 


[17] 


tha, tan /? 


CMS, = 


7 TeV, 4.7 fb" 1 


2012 limits 


[18] 


CJT 

m v — cr^o 

x XT-P 


XENON100 2012 limits (224.6 x 34 kg days) 


[21] 



Table 3. Summary of experimental constraints that enter in the computation of the likelihood function. 
The upper part lists the observables for which a positive measurement exists. For these quantities mean 
values, experimental (a) and theoretical (r) uncertainties are given, which are added in quadrature in the 
Gaussian likelihood. <5a^ USY = ajj xp — a^ M corresponds to the discrepancy between the experimental value and 
the SM prediction of the anomalous magnetic moment of the muon (g — 2)^; mh stands for the mass of the 
lightest Higgs boson, for which we use the latest CMS constraint [19]. The lower part shows observables for 
which only experimental limits currently exist, including recent limits from LHC SUSY searches [17, 18], and 
constraints on the dark matter mass and spin-independent cross-section from the XENON100 direct detection 
experiment [21]. 

that would normally not have been saved in the posterior chains (as they belong to rejected 
steps in the sampling). This results in a combined total of 348M (205M) samples for the 
cMSSM (NUHM) scans including the <5a^ USY constraint, and 323M (147M) samples for scans 
excluding the <5a^ USY constraint, out of which the profile likelihood results are obtained. For 
the cMSSM, this is a factor of ~ 100 more than our previous works [4, 5], and a factor of 
> 3 more than the frequentist global fits analysis presented in Ref. [11] (for the NUHM our 
resolution is comparable). 

As another check of the robustness of our scanning procedure, we run 10 scans in parallel 
(for both flat and log priors) for each case we consider and we compare the resulting profile 
likelihood (and best-fit points) across the different scans, and between each scan and the 
merged samples obtained from all the scans together. We have found that while each scan 
is more noisy than the combined samples (as expected), our results are consistent across all 
the scans. 

The total computational effort for the various cases considered is approximately 22 (13) 
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CPU years for the cMSSM scans, and 72 (61) CPU years for the NUHM scans including 
(excluding) the <5a^ USY constraint, for a total computational time of approximately 168 CPU 
years. 



3 Results for the cMSSM 

We begin by showing in Fig. 1 the combined impact of all present-day experimental con- 
straints, including the ATLAS exclusion limit on the cMSSM mass parameters derived from 
5.8 fb _1 total integrated luminosity, the latest measurement of the mass of the lightest Higgs 
boson from CMS and the most recent XENON100 limit on the dark matter parameters, on 
the cMSSM parameter space. Results are shown in the (m^2,mo) plane (left), the (tan /3, Aq) 
plane (centre) and the (m^o,<7?o_ ) plane (right). The top row depicts the posterior pdf for 
flat priors, the central row the posterior pdf for log priors, and the bottom row the profile 
likelihood. In each panel, the 68%, 95% and 99% credible/confidence intervals are shown. For 
comparison, blue/empty contours show the favoured regions from Ref. [4], which included all 
experimental constraints available in December 2011, previous to the discovery of the Higgs 
boson and the latest XENON100 results. 

3.1 Impact of ATLAS exclusion limit 

The latest ATLAS exclusion limit cuts further into the low-mass regions of the cMSSM, now 
excluding the entire h-pole region which was previously viable and is visible as a vertical 
region of empty/blue contours at small mi/ 2 in the (my 2 ,mo) plane. Additionally, contours 
are further pushed towards higher values of mi/2- Focusing on the (nii/2> m o) plane, we 
observe that the posterior pdf for the log prior exhibits a bimodal nature, with two connected 
favoured regions, one corresponding to the A-funnel (AF) region at high masses, and one to 
the stau-coannihilation (SC) region, corresponding to ~ 800 — 1000 GeV and mo ~ 

300 — 400. In the SC region the lightest stau is only slightly heavier than the neutralino LSP. 
Therefore, in this region the relic density of the neutralino is reduced by neutralino-stau 
coannihilations in the early universe, in agreement with the WMAP constraint. The AF 
region is characterised by a relatively light pseudoscalar Higgs, with m^,o ~ 2?n^o, which can 
mediate resonant annihilations of the neutralino LSP, making it easier to satisfy the WMAP 
relic density constraint. 

The SC region is also where our overall best-fit point is located (see below for further 
discussion). In contrast, the SC region is not favoured in the flat prior scan. The flat 
prior gives a much larger statistical a priori weight to regions at large values of the mass 
parameters, so that the corresponding posterior pdf is strongly affected by volume effects at 
high masses, and therefore peaks there. The log prior scan explores the low mass regions in 
much more detail (see [27] for a detailed discussion), so that the posterior distribution for log 
priors also favours the SC region at small values of mo and m^2, still allowed by the LHC 
exclusion limit. 

The profile likelihood function (bottom left panel) favours the SC region around the 
best-fit point, and is much more strongly localised than the Bayesian pdf. Small scalar 
and gaugino masses are strongly favoured, values of niQ > 1 TeV are disfavoured at 99% 
confidence level. Compared to the results before inclusion of the Higgs mass measurement 
(blue/empty contours), the profile likelihood contours are confined to a much smaller region. 
This is the result of two effects: pre-Higgs constraints favoured relatively low Higgs masses, 
with a best-fit value nih = 115.6 GeV (see [4]). In the region now favoured by the profile 
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Figure 1. Constraints on the cMSSM including all available present-day data (WMAP 7- year, ATLAS 5.8 

fb _1 SUSY null search and CMS Higgs mass value, XENON100 2012 direct detection limits, and others - see 
Table 3). Black, filled contours depict the marginalised posterior pdf (top row: flat priors; middle row: log 
priors) and the profile likelihood (bottom row), showing 68%, 95% and 99% credible/confidence regions. The 
encircled black cross is the overall best-fit point, obtained from about 350M likelihood evaluations. Blue/empty 
contours show constraints as of Dec 2011, before the latest LHC and XENON100 results, for comparison (from 
[4]). In the left-hand plots, the dashed/green line shows the current LHC 95% exclusion limit, while in the 
plots on the right the red/solid line represents the 90% XENON100 exclusion limit (from Ref. [21]) rescaled 
to our fiducial astrophysical dark matter distribution. We also show the expected reach of XENON1T as a 
red/dashed line. 



likelihood the constraint on rrth can be fulfilled due to the maximal mixing scenario (see 
below). Additionally, small masses are strongly favoured by the constraint on the anomalous 
magnetic moment of the muon, which can only be satisfied in the SC region, since in the AF 
region the value of g — 2 tends to 0, producing a ~ 3a discrepancy with the data. Largely 
due to these two constraints, the profile likelihood function favours a small region at low mo, 
m^ 2 that achieves high likelihood values. The importance of the g — 2 constraint in confining 
the profile likelihood function to small masses will be discussed in detail in section 3.5. 





LHC 2012 
+ XENON100 


LHC 2012 


cMSSM 


+ XENON100 

w/o Sa s u USY 


Best-fit cMSSM parameters 


m [GeV] 


389.51 


321.08 


mi/2 [GeV] 


853.03 


839.84 


A [GeV] 


-2664.79 


-2163.28 


tan /3 


14.50 


13.48 


Best-fit nuisance parameters 


M t [GeV] 


173.539 


173.779 


m b (m b ) Ms [GeV] 


4.234 


4.202 


[a em {M z )^ s ]- 1 


127.956 


127.970 


a s (M z ) MS 


0.118 


0.119 


Pioc [GeV/cm 3 ] 


0.417 


0.373 


^isr [km/s] 


224.6 


228.8 


v esc [km/s] 


549.0 


518.2 


u d [km/s] 


276.7 


274.7 


/tm x 10 2 


2.708 


2.693 


hd x io 2 


3.814 


3.886 


fr s 


0.372 


0.421 


Best-fit observables 


m h [GeV] 


123.8 


123.3 


m^o [GeV] 
5a^ sy x 10 10 


362.7 
25.47 


355.9 
3.68 


BR(B -> V s7 ) x 10 4 


3.11 


3.17 


A _ x 10 2 


8.35 


5.34 


Bi2(5 fl -»• x 10 9 


2.94 


2.94 


£i?(D s -»• tv) x 10 2 


5.10 


5.10 


5i2(L> s -»• (iv) x 10 3 


5.24 


5.24 


5i?(L> -> pv) x 10 4 
[Pb] 


3.85 

7.0 x 10~ n 

2.1 x 10" 9 


3.85 
1.1 x 10~ 10 

5.5 x 10~ 9 




0.1203 


0.1105 



Table 4. Best- fit modei parameters (top section), nuisance parameters (central section) and derived 
observables (bottom section) in the cMSSM. The column "LHC 2012+XENON100" denotes the case 
where all data, including LHC and XENON100 data, are applied in the analysis; astrophysical and 
hadronic nuisance parameters have been included in the scan and profiled over. The column "LHC 
2012 + XENON100 w/o Sa^^" is for the case where all constraints except for the Sa^ USY constraint 
are applied. 

3.2 Impact of latest Higgs mass measurement 

The new measurement of the Higgs mass has a significant impact on the cMSSM parameter 
space. For both choices of priors, large regions of parameter space previously favoured at 68% 
marginal posterior probability are ruled out by this constraint. Contours are pushed towards 
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LHC 2012 
+ XENON100 


LHC 2012 


cMSSM 


+ XENON100 

w/o baf SY 




Gaussian constraints 


SM nuisance (4 parameters) 


0.437 


1.047 


Astro nuisance (4 parameters) 


0.106 


0.725 


Hadronic nuisance (3 parameters) 


0.039 


0.237 


M w 


1.434 


1.424 


sin 2 9 e ff 


0.039 


0.039 


5a s ^ USY 


0.153 


N/A 


BR{B -> X s7 ) 


1.251 


0.933 


AM Bs 


0.133 


0.132 




1.376 
5.219 


1.377 
0.946 


BR{B u ^tv) sm 

A - 


BR(B-^Dtu) 
BR(B^Deu) 

Rm 


0.801 
0.020 


0.801 
0.020 


BR(D S -»• ri/) 


0.540 


0.541 


£i?(L> s -»• //z/) 


1.466 


1.466 


£i?(L> -> /«/) 


0.008 


0.008 




0.500 


0.001 




0.925 


1.413 




0.028 


0.030 




Exclusion limits 


XENON100 


0.633 


0.572 


LHC 


0.0 


0.0 


Sparticles (LEP) 


0.0 


0.0 


Total x 2 


15.11 


11.71 


Total Gaussian \ 2 (dof) 


14.48 (11) 


11.14 (10) 


Gaussian x 2 /dof 


1.32 


1.11 


p-value (Gaussian constraints only) 


0.21 


0.35 



Table 5. Breakdown of the total best-fit x 2 by observable for the cMSSM, for both the case where all 
data, including LHC 2012 constraints and XENON100 data, were applied, and nuisance parameters 
were profiled over (left column), and when all data except the Sa^ USY constraint were applied (right 
column) . 

larger values of m^2, the posterior pdf with flat priors excludes values of rrii/2 < 1 TeV at 
99% level. Both posterior distributions now favour much larger values of mo, with the 68% 
(95%) contour touching the prior boundary for flat (log) priors. The strong impact of the 
LHC Higgs constraint on this parameter space is expected, since in the cMSSM smaller values 
ofm/, are otherwise favoured [4]. Larger values of can be achieved by radiative corrections. 
At one- loop level is most sensitive to the stop mass, which is mainly determined by the 
value of mx/2, so that large can lead to large values of m/,, satisfying the experimental 
constraint. The value of uiq is much less important for , so that relatively low values of 
mo are still allowed by the new Higgs constraint. 
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Figure 2. Favoured regions in the cMSSM in the Xt/Ms vs rrih plane (all data included; from left to right: 
posterior pdf with fiat and log priors, and profile likelihood). The maximal mixing scenario (Xt/Ms ~ 2.44) 
is realised in the stau co-annihilation region (where the best-fit point is located, encircled black cross), while 
it can not be achieved in the A-funnel region (where Ms is larger and \Xt/Ms\ is reduced). 



A second possibility to achieve rrih ~ 125 GeV is the so-called maximal mixing scenario. 
If the stop mixing parameter Xt approaches a value ^/6M$, with Mj = 0.5(m^ + fnf 2 ), the 
contribution of the stop to rrih is maximised, and hence the Higgs mass increases (see for 
instance [54]). This is illustrated in Fig. 2, where we plot the ratio Xt/Ms vs rrih, showing 
that indeed the highest Higgs mass values are found in the maximal mixing region, where 
\X t /M s \ w VQ ~ 2.45. In the cMSSM, this effect can be achieved in the low-mass SC region, 
but is very difficult to achieve for large m±/ 2 [55]. Therefore, the posterior pdf with flat 
priors (left-hand panel of Fig. 2), which favours the AF region at larger values of m^ 2 , is 
concentrated at moderate values of Xt/Ms- A relatively large Higgs mass of rrih ~ 122 
GeV can still be achieved in this region, due to the large stop masses. Larger values of m/,, 
in agreement with the experimental constraint, would require very large values of the stop 
masses. However, the AF region can only be achieved at intermediate values of m^ 2 that 
are not large enough to result in very high stop masses leading to rrih ~ 126 GeV. 

In contrast, the posterior with log priors (central panel) and the profile likelihood func- 
tion (right-hand panel) favour much larger values of Xt/Ms, corresponding to the mode at 
small min, and a small number of fine-tuned points exist for which the maximal mixing 
scenario is realised. Nevertheless, without an additional contribution from large stop masses, 
it is difficult to achieve the measured value of the Higgs mass, and our best-fit value is 
rrih = 123.8 GeV, which is compatible with the experimental constraint (rrih = 125.8 ± 0.6 
GeV) at the ~ la level only due to the inclusion of a theoretical error of 2 GeV in the 
likelihood. 

Larger values of rrih could in principle be achieved in the Focus Point region, however, 
inside our prior range for rriQ this region is strongly disfavoured by the XENON100 constraint. 
Therefore, a value of rrih ~ 126 GeV cannot be achieved in the cMSSM within the prior ranges 
adopted for this work 2 . 

A similar pattern as in the (mi/2> Wo) plane is observed in the (tan /?, ^4o) plane (central 
panels in Fig. 1). Previously favoured regions shrink significantly due to inclusion of the 
new Higgs constraint. The posterior pdf with flat priors spans a large range of Aq values, 
with a preference for positive Aq. Large values of tan/3 are favoured, as required for the AF 

2 For an analysis of the impact of the measurement of the Higgs mass on the cMSSM at masses outside our 
prior ranges see e.g. Ref. [56]. 
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region. As mentioned above, this region however corresponds to slightly lower Higgs masses 
(nth ~ 120 — 122 GeV), and is hence disfavoured by the Higgs constraint. The posterior with 
log priors shows a bimodal shape, with the mode at low tan ft corresponding to the region 
at small m 1( / 2 in the left-hand panel. Compared to our previous profile likelihood results [4], 
we observe a strong shift of the favoured region in the (tan ft, Aq) plane, with negative Aq 
now favoured. This is a consequence of the new Higgs constraint, forcing the best-fit point 
to a region of maximal mixing. This is despite the constraint on the isospin asymmetry Ao- 
disfavouring negative Aq values [57]. 

3.3 Best-fit point 

The coordinates of the best-fit point are given in Table 4 for the input cMSSM and nuisance 
parameters, as well as for some notable derived quantities. Compared to our previous best-fit 
(in Ref. [4]) we observe an upward shift of ~ 100 — 200 GeV in the mass parameters and a 
strong shift to negative Aq (as explained above), while low tan/3 remains favoured. 

The overall best-fit x 2 , broken down in terms of the contribution of each observable, is 
given in Table 5. As can be seen, by far the largest contribution to the best-fit x 2 results 
from the isospin asymmetry Ao— . The SM prediction of this quantity is already in tension 
with the experimental measurement at ~ 2a [57]. Any positive SUSY contribution to Ao- 
will therefore further worsen the fit. The SUSY contribution is minimised at positive values 
of Aq, small tan/3 and large my 2 [57]. This preference is in tension with other constraints, 
most importantly the Higgs mass measurement, which favours strongly negative Aq. This 
leads to an additional SUSY contribution to Ao_, which further increases the contribution 
of this observable to the total x 2 ■ Other contributions to the overall best-fit x 2 are much 
smaller; in particular, the best-fit point simultaneously satisfies the constraint on the Higgs 
mass, the exclusion limit from the XENON100 direct detection experiment, the relic density 
constraint and the constraint on the anomalous magnetic moment of the muon. The value 
of BR(B S — > fj, + fi~) is also in very good agreement with the new LHCb constraint on this 
quantity. 

When evaluating the p-value for the best-fit point, we only consider contributions to 
the x 2 from Gaussian-distributed observables in the likelihood. This allows us to compute 
the (approximate) p-value analytically from the corresponding chi-squared distribution with 
the number of degrees of freedom (dof) given by the number of Gaussian data points minus 
the number of free parameters (4 cMSSM model parameter, 4 SM nuisance parameters, 4 
astrophysical nuisance parameters and 3 hadronic nuisance parameters, giving a total of 15 
parameters). We find that the p-value for all data sets combined is 0.21. Therefore, even a 
strongly constrained model such as the cMSSM is not ruled out at any meaningful confidence 
level by the latest experimental data sets. 

3.4 Implications for direct detection and future SUSY searches 

The implications for the spin-independent scattering cross-section and neutralino mass are 
displayed in the right-most column of Fig. 1. As a consequence of the new Higgs measurement, 
the favoured region in the (m^o, <r?o_ ) plane is shifted towards larger neutralino masses and 

Xi Xi v 

much lower spin-independent cross-sections, especially from the profile likelihood statistical 
perspective. The best-fit point corresponds to a very small SI scattering cross-section of 

<7?o = 7 x 10~ n pb, which is challenging to explore even with future ton-scale direct 

Xi p 

detection experiments. For comparison, the expected 90% exclusion limit from the future 
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XEN0N1T direct detection experiment is indicated on the plots in the right-hand panels. 
Therefore, the discovery of a Higgs boson with rrih ~ 126 GeV renders direct detection of 
the cMSSM more difficult. However, the posterior pdf for both choices of priors displays 
a large island of probability density around m^,o ~ 600 GeV and for SI scattering cross- 
sections °£0_ p ~ 10 -9 — 10~ 8 , corresponding to the A-funnel region in the (mi/2, mo) plane. 
This region can be fully probed by future ton-scale direct detection experiments. Therefore, 
from a Bayesian statistical perspective, significant regions of the cMSSM parameter space 
currently favoured at the 95% level will be within reach of the next generation of direct 
detection experiments. This result is obtained also from a frequentist perspective if the g — 2 
constraint is excluded from the analysis, see section 3.5 below. 

The favoured values of the spin-dependent scattering cross-section in the cMSSM are 
confined to the range 6 [10~ 9 , 10~ 6 ] pb, with the best-fit point located at cr|?_ p ~ 10~ 9 

pb, and hence outside the reach of even future multion-scale direct detection experiments 
such as DARWIN [58]. Detection prospects for neutrino telescopes are similarly pessimistic. 

We do not study in detail the impact of indirect detection experiments, such as the 
Fermi Large Area Telescope. Current limits on annihilating dark matter from observations of 
dwarf spheroidals only significantly constrain low-mass dark matter m^o < 30 GeV [59] , that 
is not realised in the cMSSM. This may change as more data become available: 10 years of 
Fermi observations of dwarf spheroidals have the capability to exclude WIMPs with masses 
m^o < 700 GeV that have a thermal annihilation cross-section< av >= 3 x 10~ 26 cm 3 /s 
(assuming 30 dwarf spheroidals [59] ) 3 . Therefore, future Fermi data will have a powerful 
impact on simple SUSY parameter spaces such as the cMSSM. 

The ID marginalised pdfs and profile likelihood for some derived quantities are shown 
in Fig. 3. The ID marginal posterior distributions for log (flat) priors are shown in red/solid 
(dash-dot /blue), the ID profile likelihood function is shown in black/dashed. 

The LHC exclusion limit pushes the sparticle masses to larger values. The difference 
in the regions of parameter space favoured by the profile likelihood and the posterior pdf 
with flat priors, as well as the bimodal behaviour of the posterior pdf with log priors, are 
clearly visible for all of the sparticle masses. Due to its confinement to small mi^, the profile 
likelihood function shows a strong preference for the smallest allowed values of the lightest 
stop and sbottom masses, with m s t op i, m s bottomi < 2 TeV. The posterior pdf stretches to 
much larger values of these quantities, up to m s t op i , ^sbottoml ~ 4 TeV. Similarly, the profile 
likelihood for the lightest chargino mass favours relatively low values of m ± ~ [500, 1000] 
GeV, while the posterior pdfs reach larger values of m ± < 1500 GeV. The range of gluino 

Al 

masses favoured by the posterior pdfs is m g i U i no ~ [2, 4] TeV, while the profile likelihood is 
confined to a relatively narrow region around m g i u i no ~ 2 TeV. The profile likelihood also 
favours a very small range of m squar k ~ 2 TeV. The posterior pdf allows a much larger range 
of average squark masses, extending to m squar k < 5 TeV. 

With the upgrade of the LHC to 14 TeV collision energy, the sensitivity to heavy SUSY 
particles will be increased significantly. With a total integrated luminosity of 300 fb" 1 , 
gluinos and squarks of the first two generations with masses around 3 — 4 TeV, and the third 
generation of squarks and charginos with masses around 800 GeV will be accessible. With the 

3 While a thermal self-annihilation cross-section can be realised in several regions of cMSSM parameter 
space, in some regions much lower values are favoured. For example, in the SC region the neutralino an- 
nihilation cross-section can be much lower than the thermal value, since the relic density is reduced by 
co-annihilations in the early universe. 
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Figure 3. ID marginal pdf for flat priors (dash-dot/blue), log priors (thick solid/red) and ID profile 
likelihood (dashed/black) in the cMSSM, including all up-to-date experimental constraints. Top row, from 
left to right: neutralino mass, lightest chargino mass, lightest stop and sbottom masses. Central row: gluino 
mass, average squark mass, lightest Higgs boson mass, spin-independent scattering cross-section. Bottom row: 
BR(B S — t fi + /j,~) branching ratio, anomalous magnetic moment of the muon, dark matter relic abundance 
and isospin asymmetry. The best-fit point is indicated by the encircled black cross. 



increase of the integrated luminosity to 3000 fb" 1 , as planned with the High-Luminosity (HL- 
LHC) upgrade, the sensitivity will further improve by a few hundred GeVs [60]. Therefore, 
detection prospects of the cMSSM at the LHC remain very promising. The majority of the 
parameter space favoured by the Bayesian posterior will be accessible to the LHC with 14 TeV 
collision energy. From the profile likelihood statistical perspective an important fraction of 
the favoured gluino and squark masses will already be explored with 20 fb" 1 of data collected 
at a collision energy of 8 TeV, the remaining currently favoured region will be explored by 
the LHC operating at 14 TeV collision energy. 

The ID posterior distributions for rrih for both choices of priors peak at relatively low 
rrih ~ 121 — 122 GeV, while the profile likelihood function favours slightly larger values rrih ~ 
124 GeV. The discrepancy between the posterior pdfs and the profile likelihood illustrates 
the difficulty of satisfying the experimental constraint on rrih in the cMSSM. The posterior 
pdf takes into account volume effects, and therefore peaks at lower rrih that are easier to 
achieve. The profile likelihood is dominated by a relatively small number of points of high 
likelihood that achieve a value of rrih closer to the experimental constraint due to maximal 
stop mixing. Both distributions are offset from the measured mass rrih = 125.8 GeV, this 
value is basically not achieved in the cMSSM. 
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The ID profile likelihood for BR(B S — > fi + n~) strongly favours a small range of values 
around BR(B S — > fi + fi~) ~ 3.0 x 10~ 9 . The posterior pdf spreads over a much larger range 
of BR(B S — > fi + fj,~) > 3.0 x 10 -9 , but also peaks at relatively small values and falls of at 
larger BR(B S — > /_t + / u~). As can be seen, the new LHCb constraint on BR(B S — > fj> + fi~) has 
a very limited impact on our results. The experimental value BR(B S — > fi + /j,~) = 3.2 x 10 -9 
agrees well with the peak of both the ID posterior pdf and profile likelihood function for 
this quantity. Significantly smaller values, that would be discrepant with the experimental 
measurement, are not realised in the cMSSM. Larger values of BR(B S — > fi + ^~) can be 
achieved, and are disfavoured by this constraint. However, the previous upper limit was 
even slightly more constraining at large BR(B S — > than the current constraint. A 

more precise measurement of BR(B S — > is needed for this constraint to have a strong 

impact on the cMSSM parameter space. 

The ID profile likelihood and posterior pdf for Sa^ u are also shown in Fig. 3. Results 
for the two statistical perspectives differ strongly. The profile likelihood function peaks at 
relatively large values of 5a^ USY , in good agreement with the experimental constraint. In 
contrast, the posterior pdf for both choices of priors favours a SM-like value of the anomalous 
magnetic moment of the muon, and therefore peaks at significantly smaller values. While 
these values are in strong disagreement with the experimental constraint, they are much easier 
to achieve in the cMSSM, especially for high values of the mass parameters. The posterior pdf 
takes into account these volume effects, while the profile likelihood function, which peaks at 
the region of highest likelihood, favours values that reproduce the experimental measurement. 

From this discrepancy one can see that the constraint on 8a^ USY may have a very strong 
impact on our conclusions, especially on the profile likelihood results. In the following section 
we will discuss the impact of this constraint in more detail. 

The ID profile likelihood and posterior pdfs for the dark matter relic density are in 
good agreement with each other. Both quantities peak at the experimentally favoured value. 
In contrast, the distributions for the isospin asymmetry Ao- are discrepant with the experi- 
mental measurement Ao- = 3.1 x 10~ 2 by more than 2a. Smaller values of Ao- are difficult 
to achieve in the cMSSM, since the SM-like value is already strongly discrepant with the 
experimental measurement, and the vast majority of points in cMSSM parameter space lead 
to a positive contribution to Ao- . As can be seen, several points leading to smaller Ao- are 
found, but are in conflict with other constraints, which leads to a low likelihood value. As 
a result, the best-fit point is also located at large values Ao- = 8.35 x 10~ 2 , which explains 
the large contribution to the best-fit x 2 from the Ao- constraint, observed in the previous 
section. 

3.5 Impact of the Sa^ USY constraint 

The magnetic anomaly of the muon, a^ = \{g — 2)^ is an exciting and powerful test for new 
physics. At present, the experimental measurement and theoretical determinations are very 
precise, enough to either strongly constrain, or even give a positive signal of, new physics. 
However, the situation is still somewhat uncertain, essentially due to inconsistencies between 
alternative determinations of the SM hadronic contribution, specifically the contribution from 
the hadronic vacuum polarization diagram <5^a^. This contribution can be expressed in 
terms of the total hadronic cross-section e + e~ — > hadrons. When experimentally measuring 
this cross-section, one obtains a value for which is discrepant from current experimental 
measurements by more than 3a. In the past this discrepancy has widely been interpreted as a 
signal of new physics. In this case, the discrepancy should be cured by contributions from new 
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physics, in our case MSSM contributions. The immediate implication is that supersymmetric 
particles should have relatively small masses, in order to produce a large enough contribution, 
<5 SUSY a /i , to reconcile theory and experiment. Hence, SUSY should be found at low energies 
(accessible to LHC), mainly because of the experimental measurement of 

The lack of a SUSY signal at the LHC strongly challenges this interpretation. In previous 
works we have investigated the impact of the g — 2 constraint on the cMSSM [4, 27, 28]. We 
found that this constraint plays a dominant role in determining the favoured regions of 
parameter space in terms of the profile likelihood, and that it was the single most important 
datum disfavouring large values of mo and m^ 2 . We thus repeat the analysis presented in 
the previous section, but exclude the experimental constraint on g — 2, in order to assess the 
robustness of our conclusions with respect to omission of this contraint. 

Results of this analysis are shown in Fig. 4 (black/filled contours), where they are 
compared with the findings presented above (blue/empty contours, corresponding to the filled 
contours in Fig. 1). The posterior distributions in the (m 1 / 2 ,mo) plane are very similar to the 
results for the Bayesian analysis including the g — 2 constraint, for both log and flat priors. 
Since a large value of g — 2 can only be obtained for low gaugino and scalar masses, removing 
this constraint from the analysis leads to a slight shift of the favoured regions towards larger 
values of thq. Nevertheless, even when excluding the g — 2 constraint from the analysis, 
the posterior pdf with log priors still favours the SC region at the 68% level, although the 
posterior probability mass associated with this mode is now reduced with respect to Fig. 1. 
The smaller size of this mode is also reflected in the (tan /3, Aq) plane. The posterior pdf 
with flat priors does not find the mode at low masses and instead shows a strong preference 
for very large values of itiq and m 1( / 2 . This prior dependence of the posterior distributions 
was already observed in section 3 and is unrelated to the g — 2 constraint. Instead, this is 
a result of volume effects at high masses influencing the results for the flat prior posterior 
distribution. 

The impact of dropping the g — 2 constraint on the profile likelihood is more significant 
(bottom row in Fig. 4). Without g — 2 the A-funnel region (which is excluded at 99% level 
when the g — 2 constraint is taken into account) is now viable at the 95% confidence level. 
Large values of mo remain viable, and the 95% confidence region extends all the way to 
the 4 TeV prior boundary. However, the best-fit point and thus the region most favoured 
from the profile likelihood statistical perspective, is still found in the SC region. Indeed, the 
coordinates of the best-fit (as displayed in Table 5) are only slightly different from before. 
We also notice that the profile likelihood results agree fairly well with the posterior pdf 
obtained with the log prior. From the hypothesis testing perspective, the best-fit point has 
X 2 /dof = 1.11, corresponding to a p-value of 0.35. Therefore, when dropping the g — 2 
constraint, the cMSSM best-fit point remains perfectly viable in light of the latest data. 

The bimodal behaviour of the profile likelihood function is also apparent in the (m^o , o"|o_ p ) 
plane (right-most panels in Fig. 4) . In contrast to the results including g — 2, a large region 
corresponding to large tan (3 values and larger neutralino masses of m^o > 500 GeV is now 
allowed at 95% confidence level. Compared to the case with g — 2, values of Op ~ 10~ 8 pb 
for a neutralino mass m^o ~ 650 GeV are now within the 95% CL. In this case, the latest 
XENON100 results actively constrain this region, which is bounded from above by the direct 
detection limit. This entire AF region will be explored by the XENON1T experiment. The 
spin-dependent cross-section remains, as before, out of reach of proposed future experiments. 

The ID marginalised pdfs and profile likelihood for the same derived quantities as in 
Fig. 3 are shown in Fig. 5. The posterior pdf with flat priors is qualitatively very similar 
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Figure 4. Constraints on the cMSSM as in Fig. 1, but now excluding the g — 2 constraint (black/filled 
contours). From top to bottom: posterior pdf with flat priors, posterior pdf with log priors and profile 
likelihood. The encircled black cross is the overall best-fit point. For comparison, blue/empty contours do 
include the g — 2 constraint (as the filled contours in Fig. 1). Dropping the constraint on the anomalous 
magnetic moment of the muon hardly changes the best-fit point, but does open up the A-funnel region in the 
profile likelihood analysis. 



to before. This is to be expected, since large ttt-o and rrii^ which lead to values discrepant 
with the experimental constraint on g — 2 were already favoured in the previous analysis. 
Results for the posterior pdf with log priors also qualitatively agree well with previous results, 
although a significant shift of posterior probability from the low-mass (in the SC region) to 
the high-mass (in the AF region) mode can be observed for all sparticle masses. Results 
for the ID profile likelihood functions are, instead, very different from before. While the 
g — 2 constraint confined the profile likelihood for the sparticle masses to narrow regions 
around the best-fit point at small masses, after dropping this constraint the profile likelihood 
function is much more spread out (although the best-fit point hardly changes). While the 
increased preference for larger sparticle masses for both the posterior distributions and the 
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Figure 5. As in Fig. 3, ID marginal pdf for flat priors (dash-dot/blue), log priors (thick solid/red) and 
ID profile likelihood (dashed/black) in the cMSSM including all current experimental constraints except the 
g — 2 constraint, for some derived quantities of interest. The best-fit point is indicated by the encircled black 



profile likelihood function slightly worsens detection prospects in this case, the majority of 
the favoured regions are still accessible for the LHC operating at 14 TeV collision energy. 

The posterior pdfs for are almost identical to the distributions shown in Fig. 3 (up 
to numerical noise). The profile likelihood function is shifted to slightly larger values of rrih, 
favouring values of rrih = [123, 126] GeV. This shows that when dropping the g — 2 constraint 
from the analysis it is much easier to achieve large Xt, fine-tuned to lead to values of 
that are in reasonably good agreement with the experimental constraint. 

The posterior pdfs for BR(B S — > /i + fi~) are qualitatively very similar to the results in 
the previous section. The profile likelihood function however now spreads over a much larger 
range of values and closely resembles the shape of the posterior pdf with log priors. 

The posterior distributions for 5a^ USY are almost identical to the distributions shown 
in Fig. 3. The posterior pdf with log priors still displays a bimodal shape, however now the 
majority of the posterior mass is found in the mode at lower values of 5a^ USY . Results for 
the ID profile likelihood function change significantly. In the absence of the experimental 
constraint on <5a^ 5y the profile likelihood is no longer pushed to large values. Instead, it 
now spreads over a large range of g — 2 values and is in much better agreement with the values 
favoured by the posterior pdfs (although it still peaks at slightly larger values of 5a^ USY ). 

Both the posterior pdfs and the profile likelihood for the dark matter relic density are 
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in good agreement with the experimental constraint, similarly to the results of the analysis 
including the g — 2 constraint, although now the profile likelihood appears less Gaussian. 
Likewise, the posterior pdfs for the isospin asymmetry, shown in the left-hand panel of the 
bottom row, are in good agreement with the previous results. However, the ID profile 
likelihood for Ao- shows a very different behaviour. This distribution is dominated by a 
small number of strongly fine-tuned points in parameter space that achieve a negative SUSY 
contribution to Ao_, so that this quantity is in better agreement with the experimental 
constraint, while also reproducing other measurements. These points show up as 'spikes' 
of high likelihood in all panels of Fig. 5. In the cMSSM, a large amount of fine-tuning is 
required to satisfy the g — 2 constraint. In the absence of this constraint there is significantly 
more freedom to find regions in parameter space that are strongly fine-tuned to satisfy other 
experimental constraints, such as Ao-. The presence of a small number of fine-tuned points 
achieving a very high likelihood value also explains the small size of the 2D 68% confidence 
level in Fig. 4. The 95% and 99% C.L. also receive contributions from points at higher Ao-, 
and are therefore much more spread out. 

3.6 Comparison with other analyses 

This is the first study that includes the new LHCb measurement of BR{B S — > the 
most up-to-date CMS constraint on the mass of the Higgs boson (derived from a combination 
of 5.1 fb" 1 = 7 TeV data and 12.2 ftr 1 y/s = 8 TeV data), and the 5.8 fb - 1 integrated 
luminosity exclusions limit from ATLAS SUSY searches. 

Other recent global fits analyses of the cMSSM can be found in Ref. [11-13]. In Ref. [11] 
a frequentist analysis of the cMSSM is presented that can be compared to our profile likelihood 
analysis. In Ref. [11] the best-fit point is found in the SC region, in agreement with our 
findings. However, the p- value is significantly lower than ours (0.085 compared to 0.21), 
perhaps as a consequence of the lower resolution of their scan, which could not find a better 
best-fit. As a consequence, the profile likelihood contours in Ref. [11] encompass the AF 
region at ~ la level even when including the g — 2 constraint. A Bayesian analysis of the 
cMSSM is found in Ref. [12]. The Bayesian posterior results agree qualitatively with ours, 
although there are important quantitative differences due to our more constraining data sets. 
Also, we found the Focus Point region to be disfavoured at > 99% level by XENON100 data, 
which are not included in Ref. [12]. In contrast to this work, Ref. [12] presents results for only 
one choice of priors (log priors), and does not discuss the prior dependence of the results. 
Their best-fit point is found in the AF region, which is excluded at 99% C.L in our profile 
likelihood analysis. We notice that the MultiNest settings used in Ref. [12] are inadequate 
to achieve a reliable exploration of the profile likelihood, as demonstrated by Ref. [25], which 
means that the best-fit point found in [12] is unlikely to be reliable. Global fits of the 
cMSSM from the Fittino group can be found in Ref. [13], for both the Bayesian and the profile 
likelihood statistical perspective. This analysis does not include the experimental constraint 
on the mass of the Higgs boson from the CMS or ATLAS collaborations. A discussion of a 
potential Higgs discovery at mh ~ 126 GeV is provided, that qualitatively agrees with our 
results. However, since the experimental and theoretical errors on assumed in Ref. [13] 
are significantly larger than in this work, and the most recent limit from the XENON100 
experiment is not included in the analysis, our results are not directly comparable. 
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4 Results for the NUHM 

4.1 Impact of all present-day experimental constraints 

The constraints on the NUHM parameters obtained from all present-day data sets, including 
the LHC 5.8 fb -1 exclusion limit and the constraint on the mass of the lightest Higgs boson, 
are shown in Fig. 6. This figure also compares the results to the constraints obtained without 
inclusion of the XENON100 2012 results. Results are shown in the (w,i/2 ) ?71 o) plane, the 
(tan/3, Aq) plane and the (m^/x) plane, in terms of the posterior pdf for flat (log) priors 
(top and central row) and of the profile likelihood (bottom row). 

The posterior distributions for both log and flat priors strongly favour large values of 
the mass parameters, especially mi/2- The LHC exclusion limit has essentially no impact on 
the parameter space, since values of m^ 2 above the LHC limits are already favoured by other 
constraints. The rrih constraint plays a dominant role. As in the cMSSM, the mass of the 
lightest Higgs boson scales with mu2- As a result, values of ?n 1 / 2 < 2 TeV are excluded at 
99% level. Since the mass of the lightest Higgs is not very sensitive to mo, almost the entire 
prior range of mo is allowed, with the exception of very small values. Posterior distributions 
for the two different prior choices agree quite well, although contours for the posterior pdf 
with flat priors are shifted towards larger values of mo, due to volume effects. 

In the (tan /3, ^4o) plane (top central panels) only very limited constraints are placed 
on the parameters. The 68% credible interval spans almost the entire prior range of Aq, for 
both the posterior pdfs with log and flat priors. Values of tan /3 < 40 are favoured, but larger 
values tan /3 ~ 50 are still allowed at 99% level. 

The Higgsino mass parameter [i is strongly constrained to values [i ~ 1 TeV, indepen- 
dent of the choice of prior. As the neutralino is Higgsino-like (as shown below), the WMAP 
relic abundance fixes its mass and hence (i. In contrast, m^ is almost unconstrained within 
the prior range, only very small values of rriA are ruled out, mainly as a consequence of the 
BR(B S — > constraint [61]. For the flat prior, the 68% contour stretches to larger 

values of uia, as expected, since the flat prior gives a large a priori statistical weight to large 
masses. The log scan shows a slight preference for values rriA < 3 TeV, but the 95% region 
still touches the upper prior boundary of tua- 

Turning now to the profile likelihood (bottom row of Fig. 6), results are qualitatively 
similar to the Bayesian pdf, but the favoured region in the (mi/2,w*o) plane is much more 
localised. A diagonal region at large m 1( / 2 > 2 TeV and intermediate values of mo = [1,2] 
TeV is favoured. In addition, a small island at 68% confidence survives at larger mo- Results 
in the (tan /3, Aq) and (rriA, /**) planes are also qualitatively similar to what has been discussed 
above for the posterior pdf, although slightly more localised. 

The main impact of the new XENON100 limit is to push contours towards larger values 
of mo- The regions ruled out correspond to the Focus Point region, which leads to large 
SI cross-sections and is therefore disfavoured by the XENON100 limit, as in the cMSSM 
(see Ref. [5] for a detailed discussion). The other parameters shown in Fig. 6 are relatively 
insensitive to this limit. 

The coordinates of the best-fit point are given in the left column of Table 6, the cor- 
responding contribution to the total x 2 by the individual observables is given in Table 7. 
The best-fit point is found at large m^ and intermediate mo- It corresponds to a slightly 
negative Aq, and a small tan/3 value. However, given the large extent of the 68% confidence 
region in the (^4o>tan/3) plane, there are many other values of Aq (and up to tan/3 < 40) 
that deliver a comparably good quality of fit. As was already the case for the cMSSM, the 
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Figure 6. Constraints on the NUHM parameters including all available present-day data (WMAP 7- 
year, LHC 5.8 fb" 1 SUSY null search and Higgs detection, XENON100 2012 direct detection limits). Black, 
filled contours depict the marginalised posterior pdf (top row: flat priors; middle row: log priors) and the 
profile likelihood (bottom row), showing 68%, 95% and 99% credible/confidence regions. The encircled black 
cross is the overall best-fit point, obtained from approximately 200M likelihood evaluations. Blue/empty 
contours show constraints without the latest XENON100 results, for comparison. In the left-hand plots, the 
dashed/green line shows the current LHC 95% exclusion limit. 



largest contribution to the overall best-fit x 2 results from the isospin asymmetry Ao_ (see 
below). Other experimental constraints are in good agreement with the best-fit point. The 
p-value corresponding to this best-fit point is 0.26, so that from the frequentist statistical 
perspective this model cannot be ruled out. 

4.2 Implications for direct detection and future SUSY searches 

The favoured regions in the (m^,o,a-\_ ) plane are shown in Fig. 7. Before inclusion of 
the new XENON100 limit, the posterior distributions (blue/empty contours) favour spin- 
independent cross-sections in the range cr?J_ = [1CU 7 , 10~ 9 ] pb, and a relatively small range 
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Figure 7. Constraints on the neutralino mass and SI scattering cross-section in the NUHM, including all 
available present-day data (WMAP 7-year, LHC 5.8 fb" 1 SUSY null search and Higgs detection, XENONIOO 
2012 direct detection limits). Black, filled contours depict the marginalised posterior pdf (left: flat priors; 
middle: log priors) and the profile likelihood (right), showing 68%, 95% and 99% credible/confidence regions. 
The encircled black cross is the overall best-fit point. Blue/empty contours show constraints without the 
latest XENONIOO results. The 90% XENONIOO exclusion limit (from Ref. [21]) is shown as a red/solid line. 
Ref. [21] only shows the limit for m^o < 1000 GeV; we show the extension of this limit to higher WIMP 

masses as a red/dashed line. We also show the expected reach of XENON1T as a red/dashed line. 



of neutralino masses around m^o ~ 1 TeV, as a consequence of the Higgsino-like character of 
the neutralino. As can be seen by comparing the blue and the black contours, the XENONIOO 
2012 limit (red/solid line) rules out part of this otherwise unconstrained region. The picture 
is similar in terms of the profile likelihood. We also display the expected 90% exclusion limit 
from the future XENON1T direct detection experiment. XENON1T will probe the entire 
currently favoured NUHM parameter space, independently of the statistical perspective. 
Therefore, direct detection prospects for this model remain very good given all present-day 
experimental constraints. 

As in the cMSSM, the spin-dependent cross-section remains out of reach even for future 
multi-ton scale detectors. The favoured region spans the interval cr??_ € [10~ 5 - 5 , 10" 6 - 5 ] pb, 
with the best-fit point found at the bottom end of the range. 

In Fig. 8, we show the ID posterior pdf and profile likelihood for some observables and 
derived quantities. The first six panels show some sparticle masses of interest. Both the 
lightest neutralino and lightest chargino masses are highly concentrated around 1 TeV. The 
favoured masses of the lightest stop and bottom are m st0 pi rj 5000 GeV and m s b ttomi ~ 
6000 GeV, respectively. The favoured gluino and average squark masses are even larger, 
"Vuino) ^squark ~ 7000 GeV. The favoured sparticle masses are far beyond the current reach 
of the LHC, and will not be accessible to the LHC operating at 14 TeV collision energy, 
nor the HL-LHC upgrade. This is true for both the posterior pdf and the profile likelihood 
function, which are in excellent agreement. Therefore, detection prospects of the NUHM at 
colliders are dim, and for discovery of this model alternative search strategies, such as direct 
detection experiments, have to be relied on. 

As can be seen from the ID distributions for m/,, in the NUHM a Higgs mass rrih ~ 126 
GeV can easily be realised, and both the Bayesian pdfs and the profile likelihood function 
peak at the experimentally measured value. The reason why rrih ~ 126 GeV is easily achieved, 
while this value is disfavoured in the cMSSM, is that in the NUHM much larger values of 
mi/2 are allowed, leading to larger stop masses, and thus larger values of rrih- The favoured 
regions in NUHM parameter space in the (Xt/Ms,rrih) plane are shown in Fig. 9. As can 
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NUHM 


LHC 2012 + 
XENON100 


LHC 2012 + 
XENON100 

w/o 5al USY 


Best-fit NUHM parameters 


m [GeV] 


1524.76 


3411.36 


mi/2 [GeV] 


3836.97 


3911.16 


A [GeV] 


-478.54 


-3519.45 


H [GeV] 


1149.27 


1132.91 


m A [GeV] 


773.47 


681.35 


tan j3 


15.37 


9.38 


Best-fit nuisance parameters 


M t [GeV] 


173.809 


173.380 


m b {m b ) Ms [GeV] 


4.240 


4.219 




127.949 


127.956 


a s (M z ) MS 


0.118 


0.117 


pioc [GeV/cm 3 ] 


0.387 


0.390 


^lsr [km/s] 


227.0 


226.5 


u esc [km/s] 


544.7 


539.5 


u rf [km/s] 


261.2 


278.9 


/t« x 10 2 


2.788 


2.819 


frd x 10 2 


3.699 


3.985 




0.360 


0.359 


Best-fit observables 


m h [GeV] 


126.1 


126.2 


m^o [GeV] 
Sa^ SY x 10 10 


1169.1 
27.92 


1153.0 
0.19 


5^(5 -> X^) x 10 4 


3.59 


3.64 


A _ x 10 2 


7.45 


7.40 


5i?(i? s -> / u+^~) x 10 9 


2.73 


2.71 


BJ?(.D a -> rz/) x 10 2 


5.10 


5.10 


BR(D S -> /xz/) x 10 3 


5.23 


5.24 


BR(D ->■ /ii/) x 10 4 

CT?q [pb] 

xV-p L J 
cr S ? [pb] 


3.85 
4.3 x 10" 9 

3.0 x 10~ 7 


3.85 
3.9 x 10" 9 

2.6 x 10~ 7 


^o/l 2 


0.1159 


0.1123 



Table 6. Best- fit model parameters (top section), nuisance parameters (central section) and derived 
observables (bottom section) in the NUHM. Data included in each column is as in Table 4. 

be seen, moderate values of Xt/M$ are favoured, and the maximal mixing scenario is not 
realised. By a combination of large stop masses and moderate stop mixing the Higgs mass is 
increased, so that ~ 126 GeV can comfortably be achieved. Since the cMSSM is obtained 
from the NUHM by imposing extra universality conditions, in principle the maximal mixing 
scenario could be realised at low NUHM masses, similar to what was found for the cMSSM 
in section 3.2. However, this scenario requires large fine-tuning and, due to the small stop 
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NUHM 


LHC 2012 + 
XENON100 


LHC 2012 + 
XENON100 

w/o 5a^ USY 




Gaussian constraints 


SM nuisance (4 parameters) 


0.831 


0.199 


Astro nuisance (4 parameters) 


0.343 


0.051 


Hadronic nuisance (3 parameters) 


0.215 


0.118 


M w 


1.539 


1.745 


sin 2 e ff 


0.022 


0.009 


Sa s ^ USY 


0.009 


N/A 


BR(B -> X s j) 


0.008 


0.053 


AM B , 


0.132 


0.132 


ID XL y ID ii ?T V } 

BR(B u ^tv)sm 

A - 


1.453 
3.570 


1.405 
3.500 


BR(B->Dtv) 
BR(B^Deu) 

Rm 


0.816 
0.017 


0.807 
0.019 


BR(D S -> ri/) 


0.544 


0.542 


Bi?(D s -> fJLv) 


1.472 


1.468 


BR(D -> 


0.008 


0.008 




0.141 


0.011 




0.021 


0.032 


Bi2(B s -»• ^ + /u") 


0.096 


0.105 




Exclusion limits 


XENON100 


0.070 


0.033 


LHC 


0.0 


0.0 


Sparticles (LEP) 


0.0 


0.0 


Total x 2 


11.31 


10.24 


Total Gaussian \ 2 (dof) 


11.24 (9) 


10.21 (8) 


Gaussian x 2 /dof 


1.25 


1.28 


p-value (Gaussian constraints only) 


0.26 


0.25 



Table 7. Breakdown of the total best-fit % 2 by observable for the NUHM. Data included in each 
column is as in Table 5. 

masses in this region, can only achieve values of slightly lower than the experimental 
constraint. Therefore, this region is disfavoured with respect to the high-mass region, in 
which rrih ~ 126 GeV can easily be achieved. 

The ID distributions for the gaugino fraction gj are shown in the right-hand panel in the 
central row of Fig. 8. One of the main new features in the NUHM compared to the cMSSM 
is the possibility of dark matter with a large Higgsino fraction (i.e. a small gaugino fraction 
gf < 0.3). Higgsino-like dark matter arises from the fact that in the NUHM \i is a free 
parameter, so that it can be adjusted to give the correct dark matter relic density required 
by the WMAP constraint. As can be seen in Fig. 8, current experimental constraints rule out 
the possibility of gaugino-like dark matter (gf 3> 0.5) and favour regions of parameter space 
that correspond to a strongly Higgsino-like LSP with gf < 0.1 at 99% level. For Higgsino-like 
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Figure 8. ID marginal pdf for flat priors (dash-dot/blue), log priors (thick solid/red) and ID profile 
likelihood (dashed/black) in the NUHM, including all current experimental constraints. Top row, from left 
to right: neutralino mass, lightest chargino, stop and sbottom masses. Central row: gluino mass, average 
squark mass, lightest Higgs mass boson mass, gaugino fraction. Bottom row: BR(B S —¥ fj> + [i~) branching 
ratio, anomalous magnetic moment of the muon, dark matter relic abundance and isospin asymmetry. The 
best-fit point is indicated by the encircled black cross. 
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Figure 9. Favoured regions in the (mi, X t /Ms) plane in the NUHM, including all available present-day 
data. Black, filled contours depict the marginalised posterior pdf (left: flat priors; middle: log priors) and the 
profile likelihood (right), showing 68%, 95% and 99% credible/confidence regions. The encircled black cross 
is the overall best-fit point. Blue/empty contours show constraints without the latest XENON100 results. 
Compared to the cMSSM, the maximal mixing scenario is not realised. 
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m^o, and a relatively large fi is required to fulfil the relic density constraint. 
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A mass of m x ~ 0(100) GeV underproduces dark matter in the universe, but higher masses 
m^o ~ 1 TeV can correctly reproduce the WMAP results. Increasing fx further leads to a 
larger dark matter relic density, in conflict with the experimental constraint. This explains 
the concentration of m^o around 1 TeV. 

The ID distributions for BR(B S — > are confined to a relatively small range 

BR(B S — > [i + fi~) ~ (2.5,4) x 10~ 9 . This ranges is comfortably within the la error range 
of the new LHCb measurement of BR(B S — > so that this constraint has a minimal 

impact on the NUHM parameter space. 

While for almost all other quantities the profile likelihood results are in good agreement 
with the Bayesian posterior pdfs, this agreement breaks down for the anomalous magnetic 
moment of the muon (see central left panel in the bottom row of Fig. 8). The profile likelihood 
function peaks in a region that satisfies the g — 2 constraint, while the bulk of the ID posterior 
distributions (for both choices of priors) has a strong preference for SM-like values of g — 2. 
This is mainly a consequence of the Higgs mass measurement and other constraints preferring 
large values of the mass parameters, while the g — 2 constraint is most easily satisfied at low 
masses. In the high-mass regions favoured in the NUHM the vast majority of points lead to 
a SM-like value of g — 2 so that the posterior pdfs strongly favour these values. The best-fit, 
and hence the peak of the profile likelihood, is found in a region of parameter space where the 
g — 2 constraint and the other constraints are simultaneously satisfied. This requires strong 
fine-tuning, so that only a very small number of such points are found by our scans. This 
explains the relatively lower resolution observed for the profile likelihood function in Fig. 6. 

The ID distributions for ft x h 2 agree well with the WMAP constraint. In contrast, sim- 
ilar to what was observed for the cMSSM, values of Ao- in agreement with the experimental 
measurement are difficult to realise in the NUHM, as can also be seen by the sizeable contribu- 
tion of this observable to the best-fit x 2 value in Table 7. Large negative SUSY contributions 
to Ao- are difficult to achieve in this model, so that SM-like values Ao_ ~ 8 x 10~ 2 are 
favoured. The SUSY contribution to Ao- is minimised, since the regions favoured in NUHM 
parameter space correspond to large my 2 , relatively small tan/3, and there is a preference 
for vanishing or positive Aq. 

4.3 Comparison with previous analyses 

Our results can be contrasted with an earlier study of the NUHM in Ref. [23]. In this 
study it was found that regions in NUHM parameter space corresponding to neutralino dark 
matter with a large gaugino fraction gf > 0.7, found at low values of mo and m^j, were 
strongly favoured, while regions corresponding to Higgsino-like neutralino dark matter were 
less favoured. In contrast, our results show a strong preference for Higgsino-like dark matter, 
corresponding to large scalar and gaugino masses, while gaugino-like dark matter is ruled 
out at 99% level from both the Bayesian and the profile likelihood statistical perspective. 

While a large fraction of the NUHM parameter space favoured in Ref. [23] is ruled 
out by the LHC exclusion limit, the constraint on the Higgs mass also plays a dominant 
role in our findings. As can be seen in Fig. 2 of Ref. [23], small values of rah ~ 116 GeV, 
now ruled out by the constraint, were previously favoured. In order to achieve Higgsino- 
like (instead of the previously favoured gaugino-like) neutralino dark matter one requires 
|/i| < Mi, where M\ is the soft mass of the bino. The bino mass scales with the gaugino 
mass as M\ ~ 0.4m]y 2 - Since fi ~ 1000 GeV is required by the relic density constraint, this 
explains the strong preference for large mi/2 > 2 TeV shown in Fig. 6. As mentioned in the 
previous section, large values of m^ lead to an increase of m^, making it easier to satisfy 
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Figure 10. Profile likelihood (black/filled, 68%, 95% and 99% CL) in the NUHM parameter space including 
the LHC Higgs discovery and 5.8 fb _1 exclusion limit, the latest XENONIOO limit and all other data, except 
the 8af t us "' constraint . The bottom-right plot shows the profile likelihood in the spin-independent scattering 
cross-section vs neutralino mass plane. Blue/empty contours include the Saf l USY constraint (as in Fig. 6), for 
comparison. The posterior pdfs are identical to Fig. 6, up to numerical noise, and hence are not shown here. 
The encircled black cross shows the overall best-fit point. 



the LHC Higgs constraint. Therefore, the constraint on mh strongly favours Higgsino-like 
dark matter over gaugino-like dark matter, which is predominantly found at small m\/2 (see 
Ref. [23]). 

Further differences with the study in Ref. [23] concern the statistical setup, as in that 
work the soft masses mn u , m H d of the two Higgs doublets are used as input parameters, 
rather than fj,, wia as in this work. Given the non-linear relationship between the two sets 
of quantities, flat priors in one set do not correspond to flat priors in the other, hence the 
Bayesian posterior is affected by the Jacobian of the transformation. The profile likelihood 
is, in principle, prior independent. However, the present work makes use of a much higher 
resolution to ensure a reliable mapping of the profile likelihood function than the scanning 
methodology adopted in [23]. 

4.4 Impact of the Sa^ 1 constraint 

Given the strong impact of the g — 2 constraint on the profile likelihood results it is of 
interest to compare Fig. 6 to the results obtained when excluding the g — 2 constraint from 
the analysis. We carried out a second set of scans including exactly the same constraints as 
before, except for the g — 2 constraint. The resulting 2D and ID distributions for g — 2 are 
shown in Figs. 10 and 11, respectively. The posterior pdfs are identical to what was found 
in Fig. 8, up to numerical noise (as can be verified by comparing Fig. 11 with Fig. 8), hence 
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we do not display the 2D posterior pdf for this case. 

Dropping the g — 2 constraint has almost no impact on the Bayesian results since the 
posterior pdfs including that constraint already favoured SM-like values, in tension with the 
experimental value, as discussed above. In contrast, the ID profile likelihood function differs 
strongly from the previous results. The best-fit point is shifted to much smaller values of g — 2 
(see Fig. 11), and the profile likelihood analysis now agrees much better with the Bayesian 
results: high values of the gaugino mass parameter m\ji > 2 TeV are favoured at 99% 
confidence level, while mo is almost unconstrained, with a small preference for mo > 1 TeV. 
Results in the other planes qualitatively agree with the profile likelihood analysis including 
the g — 2 constraint, but the contours are more spread out, stretching to higher values of 
tan j3 (top right panel) and spanning a larger range of itia (bottom left panel) . In the direct 
detection plane (bottom right panel), the profile likelihood contours cover a large cross- 
section range, and a significantly larger range of m^o than before. When excluding the g — 2 
constraint from the scan, results for the Bayesian posteriors for both choices of priors and 
for the profile likelihood function agree well. 

The ID profile likelihood function for the sparticle masses (see Fig. 11) are more spread 
out, but conclusions remain qualitatively similar to the analysis including the g— 2 constraint. 
Intriguingly, a second, less prominent peak in the profile likelihood is observed for values 
gf fa 1. This corresponds to the SC region, where m^o ~ m^ and the WMAP relic abundance 
is achieved via co-annihilation. In this case, the neutralino is bino-like, as in the cMSSM. 
However, this region of parameter space remains disfavoured by the Higgs mass constraint, 
even when the g — 2 constraint is dropped. Hence this secondary peak in the profile likelihood 
for gf is much lower than that providing the overall best-fit. 

Both the ID profile likelihood for BR(B S — > ^/ + /i~) and for g — 2 are now significantly 
more spread out. The profile likelihood for the relic density of dark matter is in much better 
agreement with the posterior pdfs. Since no more fine-tuning to satisfy the g — 2 is required, 
there is a lot more freedom to adjust the parameters to satisfy other constraints. For the 
same reason, a small shift to smaller values of Ao_, in better agreement with the experimental 
constraint, is observed, although qualitatively the profile likelihood for this quantity remains 
similar to the distribution obtained in the previous section. 

Taken together, these findings suggest that the g — 2 constraint strongly impacts on the 
favoured regions in the NUHM from a profile likelihood perspective, while the Bayesian pdf 
is much more robust with respect to this constraint. 

Details about the best-fit point found when excluding the g — 2 constraint from the 
analysis are given in the right-hand column of Table 6 and Table 7. The best-fit point is 
shifted towards larger mo and more negative Aq, and corresponds to a significantly lower 
value of g — 2; values of the other parameters remain similar to the best-fit point found in 
the analysis including the g — 2 constraint. The largest contribution to the total best-fit x 2 
still results from the constraint on the isospin asymmetry. The p- value of the best-fit is 0.25, 
which is almost identical to the p-value of 0.26 found when including the g — 2 constraint. 
This suggests that, while the g — 2 constraint has a strong impact on the shape of the profile 
likelihood function, the overall viability of the NUHM remains similar, independent of the 
inclusion of the g — 2 constraint in the analysis. 

4.5 Comparison with the cMSSM 

By comparing the 2D distributions for the cMSSM in Fig. 1 with the NUHM results in Fig. 6, 
phenomenological differences between the NUHM and the cMSSM become apparent. 
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Figure 11. ID marginal pdf for flat priors (dash-dot /blue), log priors (thick solid/red) and ID profile 
likelihood (dashed/black) in the NUHM, including all current experimental constraints but excluding the 
g — 2 constraint . Quantities as in Fig. 8. The best-fit point is indicated by the encircled black cross. 



In both models the measured value of the Higgs mass is difficult to achieve, so that this 
constraint plays a dominant role in determining the favoured regions of parameter space. 
However, the phenomenological consequences are very different: in the cMSSM the constraint 
on rrih leads to a strong preference for the SC region at low scalar and gaugino masses, where 
the maximal mixing scenario is realised. In contrast, in the NUHM large values of 
are favoured, since the resulting large stop masses combined with a small amount of stop 
mixing can reproduce the measured value of ra/,. Such large values of r?ii/ 2 are disfavoured 
in the cMSSM by the relic density constraint. This difference is particularly pronounced 
when comparing the profile likelihood for the two models, which favours very large values of 
m i/2 > 3 TeV in the NUHM, but is constrained to m^ 2 < 1 TeV in the cMSSM. Posterior 
and profile likelihood contours in the (tan f3, Aq) plane are much more spread out in the 
NUHM, while in the cMSSM small favoured regions can easily be identified. 

These differences also lead to very different observational consequences for the NUHM 
and the cMSSM. In the NUHM a Higgsino-like neutralino with mass m^o ~ 1 TeV is pre- 
dieted, while the cMSSM essentially always leads to gaugino-like neutralino dark matter with 
a mass of a few hundred GeV. In the cMSSM m^o ~ 1 TeV is excluded at more than 99% 
from both the Bayesian and the frequentist statistical perspective. Also, for the reasons 
explained above, in the cMSSM lower values of the SI cross-sections are favoured than in the 
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NUHM, which is especially apparent when comparing the profile likelihood functions in the 
(m f o, cx?o_, ) plane in Fig. 1 and Fig. 7. Therefore, detection prospects of the NUHM by direct 

Xi Xi p 

detection experiments are much more promising than for the cMSSM. In contrast, the LHC 
operating at 14 TeV collision energy will probe the majority of the currently favoured region 
of cMSSM parameter space, while in the NUHM very large gaugino and squark masses are 
favoured, so that most of the currently favoured NUHM parameter space will not be acces- 
sible to the LHC, even after the High-Luminosity upgrade. This difference makes it possible 
to distinguish experimentally between these two models given a positive signal at the LHC, 
or in a future direct detection experiment. 

5 Conclusions 

In this paper we have presented up-to-date global fits of two models of minimal Supersym- 
metry, the cMSSM and the NUHM, including the latest constraints from LHC SUSY and 
Higgs searches, and the XENON100 direct detection experiment. We also included for the 
first time the new LHCb measurement of BR(B S — > / u + /_t~), but found this constraint to 
have a negligible impact on both models. In contrast, the LHC constraint on the mass of the 
lightest Higgs boson has a very strong impact on the cMSSM and NUHM parameter spaces. 
Achieving large Higgs masses of m/j ~ 126 GeV in these models is difficult, and requires a 
significant amount of fine-tuning, either in the form of very heavy stops (and thus very heavy 
squarks in general), or the maximal mixing scenario. 

In the cMSSM this leads to a strong preference of the profile likelihood function for the 
stau-coannihilation (SC) region, in which the maximal mixing scenario can be realised. In 
contrast, the Bayesian posterior shows a varying degree of preference for larger values of my 2 
(in the A- funnel region), depending on the choice of priors. Conclusions about the detection 
prospects for the cMSSM depend on the statistical perspective: the posterior pdfs suggest 
encouraging detection prospects at future direct detection experiments, and reasonably good 
prospects at the LHC operating at 14 TeV collision energy. In contrast, the region favoured by 
the profile likelihood will be challenging to explore with future direct detection experiments. 
However, the profile likelihood favours very small sparticle masses, and the SC region is 
already significantly constrained by the 5.8 fb _1 ATLAS exclusion limit. This region - which 
also contains our overall best-fit point - will be further probed within the next months, and 
will fully be explored by the LHC operating at 14 TeV collision energy, leading to excellent 
prospects for either detecting or conclusively ruling out the cMSSM. A goodness-of-fit test 
using all available present-day constraints does not allow to rule out the cMSSM at any 
meaningful significance level. Although this result has to be interpreted with care, due to 
the approximation involved in adopting an asymptotic chi-squared distribution, it appears 
that the cMSSM is still not ruled out by current experimental constraints. Nevertheless, 
the region of cMSSM parameter space surviving the increasingly tight constraints set by the 
combination of all data sets is shrinking. 

In contrast to the cMSSM, where dark matter is almost exclusively gaugino-like, in the 
NUHM Higgsino-like dark matter is strongly favoured. This leads to a strong preference 
for large mi/ 2 > 3 TeV, that are disfavoured in the cMSSM due to the dark matter relic 
density constraint. As a result, very large stop masses are favoured in the NUHM that, 
in combination with a moderate amount of stop mixing, can reproduce the measured value 
of m^. The resulting preference for large sparticle masses renders detection prospects of 
the NUHM at the LHC operating at 14 TeV collision energy negative, even after the High- 
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Luminosity upgrade [60]. In contrast, the recent XENON100 limit was found to have a 
significant impact on this model, and future direct detection experiments will explore the 
entire currently favoured NUHM parameter space, from both the Bayesian and the profile 
likelihood statistical perspective. As for the cMSSM, according to the p-value of the best-fit 
point the NUHM is not ruled out at any meaningful significance by combined constraints 
from all available present-day data sets. 

The preference for Higgsino-like dark matter with m^o ~ 1 TeV at cross-sections easily 
accessible to future direct detection experiments is an important phenomenological differ- 
ence between the NUHM and the cMSSM. Detection prospects at the LHC are also very 
different for these two models, very promising for the cMSSM, while no signal at the LHC is 
expected for the NUHM. These phenomenological differences make it possible to distinguish 
between these models using future data from the LHC and direct detection experiments. This 
might require the adoption of a suitable scaling Ansatz for the local dark matter density, as 
introduced in [62]. 

While we have shown that both of these constrained models of SUSY are not ruled out 
yet in a quantitative manner (as our null results for the significance tests demonstrated), 
several highly complementary data sets now strongly impact on the cMSSM and NUHM 
parameter spaces and previously strongly favoured regions have been ruled out. The difficulty 
to simultaneously satisfy all experimental constraints in these models is becoming increasingly 
apparent, and strong degrees of fine-tuning are required to achieve satisfactory likelihood 
values. This motivates the study of more general SUSY models, such as the phenomenological 
MSSM (pMSSM) or the NMSSM, which are expected to be much more weakly constrained 
in light of current experimental data sets. The techniques used in this paper have provided a 
high-resolution mapping of the 15-dimensional cMSSM and 17-dimensional NUHM parameter 
space (including nuisance parameters), and are therefore expected to be useful in exploring 
less constrained models of SUSY with a richer phenomenology, such as e.g. the pMSSM. 
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